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Introduction 


The  dynamics  of  atomic  and  molecular  processes  at  solid 

surfaces  is  often  dominated  by  non-adiabatic  transitions  between 

quantum  states  of  the  atomic  or -molecular  system  or  between 

states  of  the  substrate.  Frequently  discussed  examples  are  those 

involving  charge  exchange  between  the  external  atom/molecule  and 

the  surface.  The  simplest  process  of  this  kind  is  the 

ionization/neutralization  of  an  atom/ion  incident  on  a  metal 
1  -7 

surface.  More  complicated  examples  are  charge  transfer  during 

8  9  1 0 

sputtering,  chemiluminescence,  electron  or  photon  stimulated 

desorption11  and  energy  transfer  via  an  intermediate  charge 

12  '3 

transfer  process.  Auger  neutralization  of  an  incident  ion* 

is  a  still  more  involved  example.  Other  important  processes  that 

may  be  described  by  curve  crossing  model  are  the  surface  induced 

de-excitation  of  electronically  or  vibrat ional ly  excited  molecule 

and,  in  some  cases,  molecular  dissociation  during  moiecuie- 

surface  collision. 

The  theoretical  treatment  of  these  processes  usually 
involves  the  classical  trajectory  approximation3-5'3-^0 
(reviewed  in  Section  II)  which  is  valid  only  when  the  kinetic 
energy  of  the  incident  atom  or  molecule  is  very  large.  For 
incident  energies  of  the  order  of  the  energy  differences  between 
relevant  quantum  levels  of  the  incident  species  or  of  the  surface 


w 


■ 
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this  approximation  fails.  While  much  ior.  scattering  work  is 
carried  out  at  high  kinetic  energy  the  low  kinetic  energy  regime 
is  very  important  in  electron  stimulated  desorption  and  in  the 
thermal  desorption  and  the  sticking  cf  alkali  atoms.  Since  a 
fully  quantum  mechanical  treatment  of  these  processes  is 
difficult  (especially  in  three  dimensions)  and  since  one  expects 
that  the  motion  of  the  incident  atom  is  essentially  classical  and 
quantum  effects  mostly  control  the  transition  between  the 
electronic  states,  one  would  like  to  seek  an  improved  trajectory 
approximation.  While  such  a  theory  may  yield  a  time  dependent 
Hamiltonian  for  the  electronic  problem  (as  does  the  trajectory 
approximation) ,  it  should  obtain  the  nuclear  trajectory  by  a 
procedure  that  will  recognize  the  continuous  occurrence  of  the 
electronic  transition  and  use  this  information  to  guide  the 
evolution  of  the  otherwise  classical  trajectory. 


In  this  paper  we  examine  the  implications  of  the  mean 
trajectory  approximation  for  charge  transfer  processes  at  metal 
surfaces.  In  this  approximation  the  "classical"  degrees  of 
freedom  (denoted  collectively  by  R)  evolve  under  the  potential  I 


* 

x.  x-  H._.(R),  where  Xj  is  the  instantaneous  amplitude 
—  1  r  j  — 

electronic  state  i  (characterizing  the  quantum  degrees 

freedom)  and  H,  .(R)  is  the  matrix  element  of  the  Kami 1 
j 

the  quantum  system  (which  depend  parametrically  on  R) 


of  the 


coni an  f 
between 


j 


i  and  j  quantum  states. 
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This  approximation  has  been  applied  by  several  workers  in 
1 5—21 

the  past  to  problems  in  molecular  dynamics.  Here  we  are 

concerned  with  its  justification  from  first  principles  as  well  as 
with  its  extension  to  situations  specific  to  surface  science. 

Thus  we  derive  the  approximation  by  two  methods:  the  first  is 
based  on  the  eikonal  approximation  while  the  second  is  obtained 
by  taking  the  extremum  of  the  action  functional  in  a  path 
integral  expression  for  the  quantum  propagator.  We  then  proceed 
to  apply  this  procedure  to  charge  exchange  processes  at  metal 
surfaces,  taking  as  a  prototype  the  ionization-neutralization 
processes  involving  an  alkali  atom  colliding  with  a  metal 
surface.  We  derive  the  equations  of  motion  which  govern  the 
•  dynamics  of  the  alkali  atom  motion  towards  and  away  from  the 
surface  as  well  as  the  transition  between  the  neutral  and  the 
ionic  state.  These  equations  which  describe  the  charge  transfer 
process  within  the  mean  potential  approximation  are  then 
supplemented  by  the  interaction  of  the  moving  atom  with  phonons 
and  electron  hole  pair  excitations  in  the  solid.  These 
interactions  give  rise  to  friction  on  the  atomic  motion  and  at 
finite  temperatures,  to  a  random  force  associated  with  the  random 
initial  populations  of  phonon  and  electron  hole  pair  states. 

The  mean  potential  approximation  is  superior  to  the 
classical  trajectory  approximation  in  that  it  responds  (in  an 
average  way)  to  the  quantum  process.  As  a  result  the  time 
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evolution  in  the  effective  trajectory  approximation  concerves  the 

average  energy,  in  contrast  to  classical  trajectory  approximation 

which  does  not.  At  the  same  time  the  equations  of  motion  in  the 

effective  potential  approximation  are  as  easy  to  integrate  as 

those  of  the  classical  trajectory  approximation.  Numerical 

results  for  alkali  ionization  and  neutralization  at  metal  surface 

within  the  effective  potential  approximation  will  be  presented  in 

22 

a  subsequent  publication. 

In  contrast  to  most  other  treatments  of  charge  exchange  at 
surfaces  (e.g.  refs.  3,  13)  we  do  not  use  the  one  electron  atomic 
orbital  picture  but  the  overall  energy  eigenstates  of  the  atom- 
metal  system.  This  enables  us  to  describe  this  process  within  a 
curve  crossing  model.  Thus  the  initial  state  for  a  neutral  atom 
colliding  with  a  metal  surface  at  zero  temperature  consists  of  a 
neutral  and  a  metal  both  in  their  ground  internal  states.  The 
state  of  the  system  following  electron  transfer  to  the  metal  is  a 
metal  with  one  extra  electron  above  the  Fermi  level  and  a 
positive  icn  in  its  ground  state.  The  energy  difference  between 
two  states  is  ®-s  -I  where  o  is  the  metal  work  function,  I  is  the 

C: 

atomic  ionization  potential  and  is  the  final  energy  of  the 
transferred  electron  measured  from  the  Fermi  level.  This  energy 
difference  becomes  smaller  at  shorter  atom-surface  distance  R 
because  I(R)  decreases  with  R.  In  the  single  electron  picture 
this  is  expressed  by  the  requirement  that  the  energy  of  the 
"ior.izaticn  level"  rises  as  R  decreases.  While  both  pictures 
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yield  similar  results  if  handled  correctly,  we  feel  that  the  use 
of  total  energy  levels  is  conceptually  preferable  to  the  one 
electron  picture  and  is  more  accurate.. 

The  paper  is  organized  as  follows.  In  Section  II  we 
describe  the  charge  transfer  process  as  a  curve  crossing.  The 
effective  trajectory  approximation  is  derived  by  using  the 
eikonal  approximations  (Section  III)  and  by  using  a  stationary 
phase  approximation  in  a  path  integral  representation  (Section 
IV) .  The  path  integral  method  permits  the  incorporation  of  the 
effects  due  to  the  kinetic  energy  operator  which  have  not  been 
derived  previously.  In  the  Sections  V  and  VI  we  show  how  to 
implement  the  resulting  equations  of  motion  for  a  charge  transfer 
problem  typical  of  surface  science,  which  involves  continuous 
manifolds  of  crossing  electronic  states,  and  phonon  or  electron- 
hole  pair  excitations. 

Finally,  in  Section  VII  we  summarize  our  results  and 
discuss  the  advantages  and  the  shortcomings  of  the  present 
treatment . 
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I I . 2  Alkali  adsorption  as  a  Charge  Transfer  Process. 

We  present  here  in  detail  the  physical  arguments  used  to 
construct  a  curve  crossing  model  for  the  alkali  adsorption  as  a 
charge  transfer  process  described  by  the  "chemical  reaction" 


Me(iD°  )  +  A  (  4» 1 
oc 


(Me')  (t(r)  +  ( <i;0  ) 


Here,  Me,  A  and  A+  denote  the  ground  state  metal,  the  neutral 

_  * 

alkali  atom  and  the  alkali  ion,  respectively.  (Me  )  denotes  the 

metal  surface  with  an  extra  electron:  the  star  indicates  that  the 

"negative  ion"  Me-  is  excited  (the  transferred  electron  is  placed 

above  the  Fermi  level) .  In  the  parentheses  accompanying  these 

symbols  we  indicate  the  orbitals  involved  in  the  charge  transfer 

process.  The  superscripts  specify  the  number  of  electrons 

occupying  each  orbital.  The  electron  to  be  transferred  is 

initially  located  in  the  atomic  orbital  and  ends  up  in  the 

emoty  metal  orbital  1  ,  whose  orbital  energy  (with  respect  to  the 

Termi  level)  is  e  .  There  is  an  infinite  number  of  final  states 

a 

_  •* 

(i.e.  an  infinite  number  of  excited,  negative  ions  (Me  )  ),  one 

for  each  oj  . 

a 


To  construct  the  energy  curves  corresponding  to  the  in:t:a- 
and  final  states  specified  by  Equation  (2.1)  we  use  the  Fermi 


snergy  or  the  surra; 


as  a  convenient  pom 


of  reference.  Sines 


in  all  the  relevant  experiments  the  sample  is  grounded. 


•n  ^  IT  n  - 
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it  is  disturbed  by  the  approaching  ion.  The  potential  surfaces 

discussed  below  correspond  to  electronic  energy  levels  of  the 

combined  metal-atom  system  and  not  to  single  electronic  orbitals. 

The  energy  difference  AV  (R)  between  the  final  and  the  initial 

ce 

states  specified  by  Equation  (2.1),  at  a  fixed  atom-surface 
distance  R  and  for  fixed  e^,  is  the  energy  I(R)  required  to 
ionize  the  atom,  minus  the  energy  4>  -  e  recovered  by  placinq 

00  CC  4  r 

the  electron  in  the  orbital  tp  .  The  ionization  potential  I(R)  is 

given  by  1  +  V.(R)  -  V  (R)  where  I  is  the  ionization  Dotential 

oo  i  n  °° 

of  the  isolated  atom  and  V^R)  and  vn(R)  are  the  ion-metal  and 

the  atom-metal  interaction  energies,  respectively.  If  the  ion 

binds  more  strongly  to  the  surface  than  the  neutral  we  have  I ( R ) 

<  I  .  In  what  follows  we  assume  this  to  be  the  case.  Note  that 
00 

the  work  function  $  appears,  rather  than  a  local  work 
2  3 

function.  This  is  because  we  define  the  ionization  potential 

I(R)  as  the  work  required  to  remove  the  electron  to  infinity; 

therefore  we  need  the  work  recovered  by  taking  the  electron  from 

infinity  to  e  . 

a 

The  energy  curves  corresponding  to  the  initial  and  the 
final  states  of  Equation  (2.1)  are  shown  in  Figure  1.  We  have 
assumed  that  the  energy  curve  corresponding  to  the  initial 
(neutral)  state  has  a  small  attractive  well  and  a  rapidly  raising 
repulsive  part.  The  ionic  curve  can  then  be  qualitatively  drawn 
by  using  the  expression 
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AV  (R)  =  I(R)-<t)  +e  =  I  +  ( V  .  ( R )  -  V  (R))-<t>  +e  ,  (2.2) 

o:  «cj  «  j.  n  »a 


which  gives  the  energy  difference  between  the  ionic  and  the 
neutral  curve.  At  large  particle-surface  distances  the 
seoaration  between  these  curves .is 


AV  (»)  =  I  -  4>  +  e 

Cc  CO  00  CX 


(2.3) 


in  what  follows  we  call  AV  («)  the  asvmDtotic  aisaatch  of  curve 

c;  *  * 

c .  Using  this  expression  we  can  rewrite  Equation  (2.2)  as 


AV  (R)  =  AV  (»)  +  (I(R)-I  ) )  . 

ot  a  Qo 


(2.4) 


'he  energy  difference  AV^( R) .  becomes  zero  at  the  distance 


given  oy 


AV 

a 

(V  - 

0  =  AV  (oo 

c: 

)  +  I(Ra)  - 

I  .  (2.5) 

00 

If  this  equation 

has  a 

solution 

the  neutral 

and  the  ionic  curv1 

cross  ( curve  I , 

Figure 

1)  at  R 

a 

and,  at  that 

distance,  the 

initial  and  the 

final  s 

tates  of 

Equation  ( 2 . 

1)  are  degenerate. 

The  charge  trass 

fer  process  resu 

lting  in  the 

creation  of  (Me 

;  ('i*)  is  most  likely  to  occur  at  Ra' 


2  4 

The  crossing  occurs  only  ir  the  "image  e:f sets"  contained 
in  the  term  I(R  )-I  are  able  to  compensate  the  asymptotic 
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s  solution  S  (R,t),  obtained  after  specifying  classical  initial 
mditions,  is  the  classical  Hamilton's  principal  function,  or 
le  classical  action. 

To  establish  the  meaning  of  the  condition?  (3.12)  and 

3.13)  we  take  matrix  elements  with  respect  to  the  amplitude  A  to 

urn  the  operator  equations  into  numerical  ones.  Since  we  expect 

:  to  behave  nearly  classically  we  use  for  the  wave  function  Eq . 

3 1 

3.10)  a  Gaussian  form  "*  defined  by 

SQ { R , t )  =  -at(R-Rt)2  +  hlnC  (3.15) 

and 

S(R,t)  =  Pt(R-Rt)  +  rt  (3.16) 

where 

Ct  =  (  2at/TTh)  1/4  (3.17 

Here  a. ,  ?_,  R  and  y  are  real  functions  of  time  and  C.  is 

t  v.  "C  v.  t 

the  normalization  constant  appearing  in  Zq .  (3.10).  The  meaning 

of  these  functions  can  be  understood  by  computing  the  following 
matrix  elements: 

(3.18) 

(3.19) 


1  (  F. ,  t )  RUi  (  R  ,  t )  dR  =  <R>=Rfa 
<?>  =  ?_ 

<!?.-?._)  2  >=  R / 4 a , 


(3.20) 
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III. 2  The  use  of  the  elkonal  method  for  generating  classical 
trajectories  in  the  one  state  problem. 

We  review  here  the  manner  in  which  the  eikonal  method 
generates  a  classical  trajectory  in  a  simple  one  dimensional,  one 
electronic  state  problem.  We  use  for  the  wave  function  the  form 


'ii(R,t)=C  A(R,t)  exp  {  iS  (  R ,  t )  /ft }  sexp  {  [S  (R,  t ) +iS  (  R  ,  t )  ] /ft  } 

(3.10) 

where  the  amplitude  A,  the  eikonal  S,  the  normalization  constant 

C,  and  the  function  S  are  real.  We  use  the  notation  R  when  the 
t  o 

position  of  the  particle  is  a  quantum  variable,  to  distinguish  it 

Inserting  '4)  (  R  ,  t )  in  the  time 
to 

P  +  ^1)  A  ,  (3.11) 

3  R 

3S/3R  is  an  unknown  operator 
representation) . 

If  we  assume  that 


from  the  classical  trajectory  R(t). 
dependent  Schrodinger  equation  leads 

[i  3SQ/3t  +  |~]A  =  H ( R , 

where  ?  is  the  momentum  operator  and 
i which  is  diagonal  in  the  coordinate 


(?  +  3S/3R)2  -  (3S/3R)2  (3.12) 

and 

3S/3t<<33/3f,  (3.13) 

c 

3  0 

the  equation  (11)  reduces  to  Kami  I  tor.- Jacobi  Equation 

3  S / 3  t  =  H(R,  33/3 R)  (3.14) 


If  tne  incident  kinetic  energy  is  small  these  two  trajectories 
are  rather  different.  Equation  (3.8)  gives  small  changes  in 
velocity,  while  Equation  (3.9)  strongly  accelerates  the  particle 
as  it  approaches  the  surface. 

The  fact  that  in  most  practical  cases  the  amplitudes  c^ft) 

and  c2(t)  are  both  non-zero  and  the  motion  of  the  nucleus  takes 

place  simultaneously  on  both  surfaces  gives  rise  to  interesting 

complications .  A  more  realistic  description  of  the  nuclear 

motion  requires  the  use  of  a  new  kind  of  force  F(c  ,c2),  which 

2 

depends  on  the  amplitudes  c^(t).  A  dependence  cn  )c.j  ,  i=l,2, 

alone  would  be  unsatisfactory  since  it  eliminates  quantum 

interference  effects.  Neither  Equation  (3.8)  not  Equation  (3.9) 

nor  any  simple  average  of  the  two  would  be  satisfactory. 

This  deficiency  of  the  trajectory  approximation  is  not 

confined  to  the  charge  transfer  problem  discussed  here;  it  is  a 

general  problem  to  be  faced  whenever  a  degree  of  freedom  that  we 

wish  to  treat  classically  is  strongly  coupled  to  a  quantum  degre 
2  9 

of  freedom  x.  A  reasonable  method  of  producing  an  mprovec 

trajectory  approximation  is  to  take  the  classical  limit  in  R 

while  creating  x  fully  quantum  mechanically.  There  are  several 

ways  of  doing  this  and  they  can  lead  to  different  "classical 

mechanics"  for  the  variable  R.  We  present  here  a  "mean 

trajectory"  approximation,  obtained  by  using  an  eikonal 

approximation  (Section  III)  and  a  path  integral  method  (Section 

IV).  A  more  elaborate  method  using  multiple  Gaussian  wave 

2  2 

packets  will  be  presented  separately. 
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=  Hij(R(t))  =  /dx<}>  .  (x,R(t)  )  He(x,R(t)  )0  (x.R(t)  ) 

(3.7) 

H  and  H22  are  the  ionic  and  the  neutral  energy  curves, 

respectively,  and  is  the  coupling  between  them.  The  matrix 

element  H^2  given  by  Equation  (3.7)  is  non-zero  because  the 

states  4>  ,  i  =  l,2  are  diabatic.  We  have  neglected  here,  as 

2  7 

customary  when  using  diabatic  states,  terms  containing  the  time 
derivatives  of  the  wave  functions  We  assume  throughout  this 

paper  that  the  dependence  of  H^fR),  i,j  =  l,2,  on  R  is  known. 


In  order  to  solve  Equation  (3.6)  we  must  propose  a  method 
for  computing  the  dependence  of  R  on  t .  If  the  incident  kinetic 
energy  is  much  higher  than  the  variation  of  Hii(R)  with  R  and  the 
difference  H22  -  H^,  we  can  use  a  straight  line  trajectory  and  a 
hard  wall  reflection  from  the  surface.  At  low  kinetic  energy, 
such  as  that  involved  in  the  processes  considered  here,  the 
choice  of  trajectory  is  rather  difficult.  If  the  particle  stays 
neutral  throughout  the  collision  process  then  R(t)  is  given  by 


mR(  t) 


3H22(R(t) ) 
3R(t) 


(3.8) 


If  the  particle  is  ionic  throughout,  then  R(t)  should  be  given  by 


mR  (  t ) 


3H1.(R(t) ) 
3R(t) 


(3.9! 
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a  classical  field.  In  what  follows  we  also  use  R  to  denote  the 
nuclear  position  appearing  as  an  argument  in  the  wave  function. 

The  state  i|i(x,R)  associated  with  the  Hamiltonian  H(x,R) 
depends  on  both  x  and  R.  If  two  electronic  states  <}>  (x,R)  and 
$  (x,R)  (2  for  the  neutral  and  1  for  the  ionic  state)  are 

sufficient  to  describe  the  system  then 

lii(x.R)  =  x1(R)4»1(x,R)  +  x2(R)4>2(x,R)  ,  (3.3) 

where  x^R)  and  x2(R)  are  nuclear  wave  functions.  In  the 

trajectory  approximation  the  total  wave  function  is  obtained  by 

renlacing  R  in  H  (x,R)  and  iJj(x,R)  with  the  trajectory  R(t)  and 
■  e 

the  nuclear  wave  functions  x^(R)  with  the  amplitudes  c^(t): 

i^(x,t)  =  c1(t)p1(x,R(t))  +  c2  ( t )  4>2  Cx,  R(  t )  )  .  (3.4) 

Inserting  Equation  (3.4)  in  the  time  dependent  Schrodinger 
equation,  and  assuming,  for  simplicity,  that 


leads  to 


Lhc  . 


=  H.  .  (R(t)  )ci(t)  +  HiJ(R(t) )Cj(t)  j , i  =  1 , 2  (3.6) 


(  i*  j 


with 


Ill  THE  MEAN  TRAJECTORY  APPROXIMATION  DERIVED  BY  USING  THE 
EIKONAL  APPROXIMATION. 

III.l  Introduction 

In  order  to  specify  in  the  simplest  manner  the  reasons  for 

the  present  work  we  first  consider  a  model  in  which  only  two 

curves  are  important.  The  dynamics  of  such  a  system  is 

considerably  simplified  by  making  the  so  called  trajectory 
3468  9-1 3  27 

approximation,  '  '  '  '  '  which  is  briefly  described  below. 

Consider  a  system  with  two  kinds  of  degrees  of  freedom  denoted  x 
and  R.  In  our  case  x  denotes  collectively  the  coordinates  of  all 
the  electrons  and  R  is  the  surface-atom  distance. 

The  full  Hamiltonian  for  this  system  can  be  written  as 

H(x.R)  =  (R2/2m)  V2  +  H  (x.R)  .  (3.1) 

e 

This  is  the  sum  of  the  kinetic  energy  of  the  nuclei  and  the 
electronic  Hamiltonian  Hg(x,R)  which  contains  the  kinetic  energy 
of  the  electrons  and  the  interaction  energy  between  electrons  and 
nuclei.  The  trajectory  approximation  uses  the  Hamiltonian 

H(x,R(t))=He(x,R(t)}  (3.2) 

where  R(t)  is  the  position  of  the  atom,  whose  time  dependence  is 
given  by  classical  mechanics.  Note  that  while  R 

appearing  in  Equation  (3.1)  is  an  operator,  the  quantity  R(t)  is 
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the  absence  of  the  approaching  atom  or  ion.  The  local  work 
function  <t>(R)  which  is  affected  by  the  approaching  species  also 
plays  an  important  role,  but  only  in  the  calculation  of  the 
matrix  elements  which  induce  the  electron  transfer.  In  this 
paper  we  assume  that  these  matrix  elements  are  known  for  all  R. 


9 


1 
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surface  residence  time  can  be  included  by  allowing  the  ion  to 
recapture  the  lost  electron  on  its  outgoing  trajectory.  A  finite 
surface  residence  time  can  be  simulated  by  adding  a  rate  of 
disappearance  of  the  electron  from  the  surface  region. 

We  conclude  this  section  with  several  remarks.  First,  even 

though  we  have  sometimes  used  a  one  electron  language,  the 

present  description  is  not  a  one  electron  model,  and  the  states 

used  here  are  correlated  electronic  states  of  the  metal-atom 

system.  Second,  even  though  we  use  the  terms  ionic  and  neutral 

throughout  the  paper,  the  model  is  not  predicated  on  the 

27 

assumption  that  the  diabatic  curves  are  purely  ionic  or  neutral 

at  all  atom  surface  distances.  Preliminary  calculations  using 

the  Generalized  Valence  Bond  and  the  Hartree-Fock  method  for  Na 

adsorption  on  Ni  clusters  show  that  the  "ionic"  states  are  not 

purely  ionic.  Nevertheless,  the  dynamic  theory  developed  here 

can  be  applied  without  modification  to  whatever  states  are 

produced  by  quantum  chemistry.  The  labels  ionic  and  neutrals  can 

still  be  used  for  those  curves  which  lead  to  alkali  ions  or 

neutrals  at  large  atom-surface  distance.  Third,  we  have  used 

throughout,  for  simplicity,  a  language  valid  for  the  one 

dimensional  case.  However,  all  our  remarks  can  be  extended  to 

three  dimensions  without  difficulties.  Finally  note  that  the 

potential  surfaces  for  the  motion  of  the  neutral  atom  and  of  the 

ion  are  characterized  by  the  local  ionization  potential  I(R)  and 

the  work  function  <b  .  The  latter  corresponds  to  the  surface  in 

00 
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(i.e.  large  e^,  large  asymptotic  mismatch  dVa(°°))  end  up  deeper 
into  the  ionic  well  than  those  occurring  at  large  R  ,  and 
have  a  smaller  chance  of  neutralization,  therefore  they  have  a 
larger  trapping  probability.  Furthermore,  the  probability  of 
trapoing  goes  up  with  4>  -I  since  this  forces  the  ionization  to 
occur  at  smaller  values  of  R  . 


This  qualitative  picture  can  also  be  used  to  understand 

coverage  effects.  As  alkali  atoms  coverage  is  increased  6  is 

'  00 

2  5 

lowered  dramatically  and  $^-1^  becomes  negative.  The  alkali 

atom  cannot  be  fully  ionized  and  adsorption  must  occur  by  usual 

chemisorption.  The  desorption  process  must  occur  exclusively 

through  neutral  desorption  until  the  coverage  is  lowered  so  much 

(by  desorption)  that  4>  -I  becomes  positive  and  both  ion  and 

neutral  desorption  become  possible.  Such  behavior  is  observed, 

2  6 

for  example  for  K  desorption  from  Ni(lll)’ 


Interesting  situations  appear  when  the  surface  residence 
time  of  the  transferred  electron  is  comparable  to  or  larger  than 
the  collision  time.  This  might  happen  either  when  the  electron 
is  transferred  into  an  empty  surface  state  with  a  lifetime  T-Tc> 
when  the  kinetic  energy  of  the  incident  particle  is  very  high,  or 
when  the  conductivity  of  the  material  is  low.  The  surf ace ’state 
can  be  included  in  the  model  presented  here  as  a  discrete  state 
with  a  finite  width  T  (T  1  is  the  lifetime  of  the  electron  in  the 


state).  The  presence  of  an  electron  with  an  infinitely  long 
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all  r>  smaller  than  rj™  given  by  4>  — n™-I  =min  &.(R). 

p  p  oo  p  oo  1 

A  newly  farmed  ion  starts  interacting  with  the  thermally 
excited  electrons  as  soon  as  it  is  produced.  The  rate  of 
neutralization  is  the  quantum  mechanical  transition  rate  to  the 
neutral  curve  labelled  by  ti  ,  multiplied  by  the  probability  that 

P 

rj  is  thermally  occupied. 

p 


We  can  now  summarize  the  overall  dynamic  picture  for  the 

case  Tr<<Tc  anc*  4>04-IOB>0 .  There  is  a  finite  probability  that  the 

incident  neutral  is  not  ionized,  in  which  case  it  will  scatter 

back  into  the  vacuum.  The  ionization  brings  about  some 

complexity.  If  it  takes  place  by  a  transition  to  the  ionic  curve 

labelled  by  a,  the  electron  is  placed  in  the  empty  metal  orbital 

from  where  it  promptly  moves  into  the  bulk  and  disappears.  If 

the  ionizing  transition  occurs  to  a  curve  ot  'for  which  A  =K-AV  <0 

oc  oc 

(i.e.  the  kinetic  energy  of  the  neutral  is  below  the  ionic  curve 
a)  the  ion  can  escape  from  the  surface  only  by  subsequent 
neutralization,  caused  by  tunneling  of  thermally  excited 
electrons  (since  the  ion  crosses  the  neutral  curves  shown  in 
Figure  2).  The  chance  of  neutralization  depends  on  the  kinetic 
energy  of  the  ion.  A  low  kinetic  energy  ion  is  deep  in  the  ionic 
well  and  can  only  undergo  transitions  on  the  neutral  curves 


located  below  it,  which  correspond  to  larger  values  of  n  , 

£ 

therefore  to  smaller  probabilities  that  n  is  thermally  occupied 

P 


by  an  electron.  Ionizing  transitions  that  take  Diace  at  small  R 

OC 


Me  (  ^  )  +  A"*"  (  (Jj  0  )  -*  N!e+{4)°)  +  A(b:)  . 

P  P 


Here  Me(iJj  )  represents  the  metal  having  a  thermally  excited 

electron  in  the  orbital  ip  ^  whose  energv  tv  is  above  the  Fermi 

level  (  t]  is  the  energy  of  the  transition  from  the  Fermi  level  tc 
p 

<|)  ),  Me~(i(j0)  is  the  metal  after  electron  removal  from  m  . 

3  3  3 


The  energy  required  for  the  process  described  by  Equation 
(2.6)  is 


AS@(R)  =  V^-K^) 


In  Figure  2  we  plot  the  ionic  curve  &, (R)  and  the  set  of 
neutral  curves  given  by 


Sn,a(R)  =  &i(R)  +  *oo"7VI(R)  • 


'We  assume,  for  simplicity  that  the  neutral  curve  is  independs 
of  R)  . 


If  r|  =0  (i.e.  we  consider  an  electron  at  the  Fermi  level) 
P 


(«) 

(« 

)=5a3~IM  and  the 

neutral  curve  is  above 

the  io 

n  z.  c 

conside 

r 

the  case  *„>!„) 

The  neutral  curve 

and 

the 

rve  &, 

do 

not  cross.  If 

b  -I  =ri  the  two  curves 

-  00  30  1  3 

coir.c 

ice 

If  1 

>0 

-I  the  centra 

1  curve  crosses  the  ioni 

c  one 

11 


(  2TT/h)  |  H12  (  R^)  |  2p  (  e^}  where  H12^Rc^  is  the  coupling  between  the 

ionic  and  the  neutral  states  at  their  crossing  distance  R  and 

a. 

p(e  )  is  the  density  of  the  ionic  states.  The  condition  r  <<t 
a  r  c 

precludes  the  neutralization  of  the  ion  by  a  recrossing  of  the 
incident  neutral  curve  on  the  outgoing  trajectory.  The  dynamics 
occurs  as  if  the  incident  neutral  curve  disappears  as  soon  as  a 
transition  to  an  ionic  curve  is  completed. 

While  the  ion  cannot  be  reneutralized  by  the  electron  which 

was  previously  transferred  to  the  surface,  neutralization  can 

occur  by  tunneling  of  thermally  excited  electrons.  The  condition 

<|>  -I  >0  precludes  the  electrons  at  or  below  the  Fermi  level 
00  00 

from  participating  in  the  neutralization  of  the  ion;  furthermore 

if  (t»oa-Ioo>>kT  only  very  few  thermal  electrons  have  enough  energy 

to  neutralize  the  ion  resonantly  and  the  neutralization  rate  is 

extremely  small.  The  behavior  of  the  ion  in  the  extreme  case  in 

which  t  < < t  and  4>  -I  >>kT,  depends  on  the  value  of  A  =  K  - 

AV^foo),  where  K  is  the  kinetic  energy  of  the  incident  atom.  If 

the  ionization  occured  by  a  transition  for  which  A  <0  the  ion 

a 

is  trapped;  if  >  0  trapping  occurs  only  if  the  ion  loses  (to 
phonons  or  electron-hole  pairs}  an  amount  of  energy  larger  than 
Aq.  Otherwise  the  ion  escapes  into  the  vacuum. 

If  4»  -  I  >0  but  is  not  much  larger  than  kT  the  ion  can 

00  00 

be  neutralized  by  tunneling  of  thermally  excited  electrons.  This 
process  can  be  described  by  the  "chemical  reaction" 
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{a),  hence  the  ionic  curves,  constitute  a  continuous  manifold,  so 
that  as  the  neutral  approaches  the  surface  the  neutral  state 
crosses  ionic  curves  over  a  large  distance  range. 

Let  us  assume  that  a  transition  to  the  ionic  state  a  has 
occured  at  R^,  while  the  neutral  was  on  its  incoming  trajectory. 
This  places  the  transferred  electron  in  the  empty  surface  orbital 
e  .  The  subsequent  dynamics  depends  on  whether  the  ion  is 
neutralized  on  the  remainder  of  the  trajectory.  Since  we  treat 
the  transition  to  the  ionic  state  cc  by  full  quantum  mechanics, 
rather  than  a  one-crossing  model,  or  perturbation  theory,  the 
transition  amplitude  contains  the  effect  of  all  the  crossings 
back  and  forth  between  the  two  states,  taking  place  while  R  is 
close  to  .  The  ion  hits  the  wall,  turns  around,  and  approaches 
again  the  point  where  is  could  be  neutralized  by  the  electron 
initially  placed  in  e^,  if  that  electron  has  not  already  moved 
into  the  bulk.  In  what  follows  we  assume  that  surface  residence 

time  t  of  the  transferred  electron  is  much  smaller  than  the 

r 

collision  time  xc<  Since,  strictly  speaking,  the  electron  is  not 
transferred  into  a  one  electron  state  but  in  a  localized  wave 
packet  centered  around  e^,  the  electron  moves  towards  the  bulk 
with  the  group  velocity  of  the  packet.  The  surface  residence 
time  of  interest  here  can  be  defined  by  tr=L/v  where  L  is  the 
displacement  required  tc  cancel  the  overlap  between  the  wave 
packet  and  the  hole  state  in  the  atom  and  v  is  the  group 

_  i 

velocity.  An  equivalent  statement  is  that  x is  roughly  given  by 
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mismatch  AV^w).  If  the  asymptotic  mismatch  is  zero  than 
Equation  (2.4)  gives  I(R  )=I  and  the  crossing  occurs  at  very 
large  values  of  R  .  As  the  mismatch  is  increased,  the  crossing 
point  moves  to  smaller  values  of  R  ,  provided  that  AVa(oo)>o. 
Obviously  if  AVa(«)<0  the  curves  cannot  cross  at  any  distance. 
Also  if  AV^oo)  becomes  too  large  the  image  field  cannot 
compensate  for  it  and  the  curves  do  not  cross.  We  denote  by 
AV*  («)  the  maximum  value  of  AVa(»)  that  can  be  compensated. 

The  ionic  curves  can  thus  be  characterized  by  any  of  the 

following  parameters:  the  asymptotic  mismatch  AV  (»),  the 

a 

crossing  distance  R^,  the  orbital  energy  e^,  or  the  index  a. 

When  AVa(«)  goes  up,  R^  gets  smaller  and  increases. 

In  discussing  the  dynamic  processes  possible  in  this  system 
of  curves  we  assume  that  we  have  a  procedure  for  generating  a 
classical  trajectory  R(t)  describing  the  motion  of  the  incident 
neutral.  This  incident  particle  can  undergo  transitions  to  any 
of  the  ionic  curves  as  long  as  they  cross  the  neutral  one.  A 
transition  to  an  ionic  state  AV  ( « )  tends  to  take  place  near  the 
crossing  point  R  .  At  any  other  point  R  there  is  an  energy 
mismatch  AV  (R),  which  means  that  after  the  transition  to  the 

OC 

curve  a  the  nucleus  must  change  suddenly  its  velocity.  Since  the 

electron  transfer  cannot  provide  a  large  momentum  transfer  such 

an  event  is  unlikely  unless  AV  (R)  is  very  small,  that  is,  unless 

oc 

R  is  very  close  to  R  .  It  should  be  kept  in  mind  that  the  states 
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We  see  that  and  are  the  expectation  values  of  the  position 
and  momentum,  y^  is  a  phase,  and  a,_  is  proportional  to  the 
quantum  fluctuation  of  the  momentum  around  its  expectation  values 
and  inversely  proportional  to  the  fluctuations  in  the  position. 


Taking  the  matrix  element  of  the  condition  (12)  with 
resDect  to  A  leads  to 


/ A ( R , t )  [P  +  3S/3R]2  A ( R , t )  dR  =  /A AdR  (3.22) 


If  we  use,  as  an  approximation  to  S  the  classical  action  Sc,  and 
the  classical  equation  P(t)  =  3S  /3R(t),  we  can  easily  perform 
the  integrals  in  Equation  (22),  to  obtain 


Hoc  << 


(3.23) 


We  have  identified  P(t)  =  3S  /3R(t)  with  the  expectation  value  of 

c 

the  momentum  P  ,  and  this  will  be  justified  later.  The  condition 
(3.12)  is  thus  equivalent  to  the  requirement  of  Eq .  (3.23)  that 

the  quantum  fluctuations  of  the  momentum  are  much  smaller  than 
its  expectation  value. 


Treating  the  condition  (3.13)  in  a  similar  manner  requires 
the  evaluation  of 


X A ( 3SQ/3t ) AdR  <<  XA ( 3  S / 3  t ) AdR 


3.24) 


Performing  the  integrals  leads  to 


j*  A  (  3Sq/3  t  )  AdR=0 


/A  ||  AdR  =  H c ( R ( t ) ,  P(t) ) 


(3.25a) 


(3.25b) 


We  have  replaced  3S/3t  with  3Sc/3t  which  in  turn  is  equal"*  i 

the  classical  Hamiltonian  H  .  Furthermore,  we  used  Eauation 

c 

(3.10)  for  S  .  Clearly  condition  (3.24)  is  alwvas  satisfied 
since  the  classical  Hamiltonian  is  positive. 


The  condition  discussed  above  establish  under  what 
circumstances  we  can  replace  the  eikonal  S(R,t)  with  the 


classical  action  Sc(R(t),  P(t; 


So  far  the  classical  mechanics 


appears  as  a  device  to  compute  the  eikonal:  the. e  is  no 
guarantee  that  the  classical  trajectory  has  any  relationship  wit 
the  motion  of  the  particle  as  given  by  the  time  dependent 
Schrodinger  equation.  For  example,  if  a„=0  the  Gaussian  wave 
packet  becomes  a  planar  wave,  and  the  condition  (3.23)  is 
fulfilled.  However,  it  is  impossible  to  describe  the  behavior  o 
a  particle  in  a  planar  wave  state  by  any  kind  of  classical 
trajectory. 


The  condition  under  which  the  trajectory  generated  by  the 
Hanti  1  ten- Jacobi  ecuation  anoroximates  the  behavior  of  the  cuar.t; 


system  can  be  obtained  by  using  the  wave  packet  defined  bv 
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Ecuat  ior.s 


.5)  ,  (3.16)  and  (3.101 


and  the  Ehrenfest  theorem 


32 


~  jV(R,t)  Ri|i(  R ,  t )  dR  = 


-;U)*(R,t)  ii)(R,t)  dR  (3.26) 
3R 


According  the  Equation  (3.18)  the  left  hand  side  of  the  above 
equation  is  dR^/dt.  Since  the  wave  function  appearing  in  the 
right  hand  side  is  localized  over  a  spatial  range  given  by 

i/2 

(2 ?»/ot  )  we  can  expand  V(R)  in  Equation  (3.26)  in  powers  of  R- 
and  this  leads  to 

s. 


cR, 


c: 


av(Rt) 


i  i3v|v 
1 


<(R-RJ 


.27) 


the  expectation  value  R  of  R  satisfies  the  classical  equation  of 
notion  if  the  secona  term  in  Equation  (3.27)  is  much  smaller  than 
the  first.  Using  Equation  (3.20)  we  can  write  this  condition  as 


!r./8a. 


3  ~V ( Rfc 


>R' 


<< 


3V(R„ ) 


3R, 


If  conditions  (3.23),  and  (3.23)  are  satisfied,  the  eikcnal 
is  given  by  the  classical  action,  and  the  trajectories  generated 
by  classical  mechanics  are  a  reasonable  description  of  the  mean 
quantum  motion  of 


the  svstem. 
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In  what  follows  we  postulate  that  the  above  procedure  can 
be  used  to  generate  a  "classical  mechanics"  for  cases  in  which 
the  variable  R  to  be  treated  classically  is  strongly  coupled  to  s 
quantum  companion  x,  and  we  apply  this  idea  to  the  case  of  a  two 
state  system. 

III. 3  The  application  of  the  eikonal  method  to  the  two 

state  problem. 

To  apply  the  method  outlined  at  III. 2  we  assume  a  wave 
function  of  the  form 


ii  (  x ,  R  ;  t )  =  A  (  X ,  R ;  t )  exp  {  iS  (  R  ,  t ) /h } 

=  (xx  (  R;  t )  ( x,  R)  +  x2  (R;  t)<t>2  (x,R) }  exp(is/h). 


(3.29) 


Here  S(R,t)  is  an  unknown  real  function  (the  eikonal)  of  the 
nuclear  coordinate  R  and  the  time  t,  and  x^tR.t)  i=l , 2  are  the 
nuclear  wave  functions  when  the  system  is  in  the  electronic  state 
d>,(x,R;  t),  i  =  l,2.  To  derive  the  appropriate  classical  limit  for 
the  nuclear  motion  we  eliminate  first  the  electronic  degrees  of 
freedom  from  the  problem  by  introducing  the  wave  function 
■1j  (  x  ,  R  ;  t )  given  by  Equation  (3.29)  in  the  time  dependent 

Schrodinger  equation,  and  by  operating  on  the  resulting  equation, 

* 

from  the  left,  with  / ax$ . ( x , R ) ,  i=l,2.  We  obtain 


||  +  ih  —]  x,  ( R ,  t )  =  (  (  1  /  2m )  {  [  ?  +  — ]  ^  H..(R))x. 

O'*.  0^1  -i.1  1 

3  R 

-  H12(R)x2 


( 3 . 30a ) 


and 


{  3 . 30b) 


[-3S/3t+ih3/3t]x2{P.;t)  =  {  (l/2m)  [P+3S/3R]“+H22(R)  }x2 

+  H2l'*)xl 

We  have  neglected  the  matrix  elements  containing  3q./3R  or 
2  2 

3  4k/3R  since  we  work  with  a  diabatic  representation  and  assume 
that  the  term  is  the  largest  coupling  between  the  electronic 

states.  The  matrix  elements  H^(R)  are  defined  by  Equation  (3.7) 
and  the  orthonormality  condition  (3.5)  has  been  assumed. 


The  classical  limit  is  taken  by  assuming  that  ih  3x^/3t 
(3S/3t)Xi  »  ( 3S/3 1 )  and  (P+3S/3R)2  ^  (3S/3R)2,  and  the 


equations  (3.30a,b)  become 


-(3S/3t)Xl  =  {(l/2m)(3S/3R}‘  +  H11>X1  +  h'12x2 


( 3 . 31a) 


~(3S/3t)x2  =  { ( l/2m) (3S/3R) ‘  +  H22}x2+H21Xl 


( 3 ,31b) 


Multiplying  Equation  (3.31a)  with  and  (3.31b)  with  x2  and 


adding  the  results  leads  to 
- 3  S / 3  t  =  ( l/2m)  ( 3S/3R) 2  +  V(R;t) 


(3.32 


wi  th 


V  (  R  ;  t )  = 


XiHi,(R)xi 


i=l  j-1 


X  j  Xj  f  (3.33) 


Equation  (3.32)  is  a  Hami 1  ton- Jacobi  equation  with  an  "effective 
potential"  V(R;t)  which  depends  on  the  potential  energies  H..(R) 


of  the  two  states,  on  the  coupling  Hiz(R)  between  them,  on  the 

*  *  _i 

probability  x.x.(Zx4x-)  that  the  particle  is  in  the  state  i, 
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(i=X,2)  and  on  the  "coherent  term"  x  j  which  introduce 

effect  of  quantum,  interference  between  the  two  nuclear  s 
As  in  Section  III. 2  classical  mechanics  appears  only  as 
to  compute  the  eikonal . 

Using  Equation  (3.32)  we  can  rewrite  the  Equations 


s  the 
fates . 
a  device 


(3.30) 


as 


ih3Xl /3t=( l/2m) (P2+  ~  -  V(R.t))x, 


3  R 


3  R  2 


*  H11(R)*1+H22(R>*2 


and 


(  3 . 34a) 


IF 


( 1 / 2m ) { P2+2  ~  ? 

3  R 


+  h 

1  3R2 


-  V  (  R  ; 


)  X. 


H22(R)X2+H21(R)*1 


( 3 . 34b! 


If  the  wave  function  is  such  that  the  momentum,  fluctuations  are 
small  compared  to  the  classical  momentum  given  by  Equation 
(3.32),  and  the  functions  of  coordinate  S(R,t),  H,.^.(R)  and  V(  R,t 
vary  smoothly  with  R  over  the  spatial  scale  set  by  the  quantum 
fluctuations  in  coordinate,  we  can  replace  P  and  R  in  Equations 
(3.34)  with  their  classical  values  ?(t)  and  R(t).  Furthermore, 
the  diagonal  matrix  element  common  to  both  equations  can  be 
eliminated  by  introducing  in  the  wave  function  'J>(x,R;t)  defined 
by  Eq.  (3.29)  the  appropriate  phase  factor.  With  these  changes 
the  (3.33)  become 
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ihi.(t)  -  Hli(R(t))xi(t)+Hi j(R(t))x.(t)  j , i=l , 2  (  i~  j  )  (3.25) 

and  the  effective  potential  is 
2  2 

V(R(t);t)  =  I  I  x*(t)x,(t)H ..(R(t))/Z  Xj  ( t ) x  j  ( t )  .  (3.25) 

i=l  j=l  J  J  i 

Here  the  nuclear  wave  function  depends  on  tine  only. 

We  can  now  summarize  the  result  obtained  by  using  the 
eikonai  method  for  the  two  state  problem.  The  quantum  amplitudes 
X^(t)  are  given  by  the  same  equations  as  those  of  the  customary 
trajectory  method  (Section  III.l).  The  "classical”  equation  of 
motion  is  however  deeply  modified  since  the  potential  energy 

given  by  Equation  (3.36)  is  neither  nor  H22,  nor  a  simple 

* 

classical  average  I  x^  x^K^tR). 


IV 


PATH  INTEGRAL  FORMULATION 


As  we  have  already  mentioned  our  strategy  in  generating 
improved  trajectory  approximations  is  to  formulate  the  problem 
quantum  mechanically  and  then  to  take  the  classical  limit  in  the 
degrees  of  freedom  whose  motion  is  to  be  described  by  a 
trajectory,  while  maintaining  a  quantum  theory  for  the  other 
variables.  The  classical  limit  can  be  taken  by  a  variety  of 
methods.  When  applied  to  all  the  degrees  of  freedom  of  a  one 
state  problem,  all  these  methods  lead  to  Newton's  equation. 
However,  the  same  methods  applied  to  a  many  state  problem  in 
which  some  degrees  of  freedom  are  treated  quantum  mechanically, 
lead  to  different  trajectory  equations.  The  acceptance  of  such 
approximations  and  the  choice  of  the  best  among  them  is  based  on 
our  prior  intuitive  expectations  of  what  such  a  mechanics  should 
be,  on  their  agreement  with  the  experiments  (where  reliable 
comparisons  can  be  made)  and  on  their  computational  advantages. 
Given  this  situation  it  seems  to  us  worthwhile  to  explore  under 
what  circumstances  various  methods  of  taking  the  classical  limit 
lead  to  the  same  results.  Here  we  show  that  within  a  path 
integral  formulation  the  theory  derived  in  Section  III  is 
recovered  if  we  use  a  coherent  state  representation  to  describe 
the  electronic  states  and  then  treat  all  degrees  of  freedom 
classically  (by  taking  the  stationary  phase  approximation). 


We  consider  the  Hamiltonian 
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N  N  N 

H  =  2_Hnti(R)tn><n|  +  Hnm(R)ln><m|  +  /  |n><n| 

n=l  n,m=l  n=l 

(4.1) 

where  | n>  are  the  electronic  states,  and  R,  P  and  M  denote 
nuclear  position  momentum  and  mass,  respectively.  This 
Hamiltonian  corresponds  to  a  diabatic  representation  in  which  the 
different  electronic  states  are  coupled  by  Hnm  and  the  coupling 
due  to  the  nuclear  kinetic  energy  operator  is  disregarded.  We 
wish  to  derive  a  path  integral  representation  for  the  propagator 

-( i/h)Ht 

K{R«K.R0«|>0I  t)  =  <R4<|e  I  R0^0>  M-2) 


where  |R4i>  denotes  a  state  of  the  system  in  which  the  electronic 
state  is  (Ji=E  xn|n>  and  the  nuclear  state  is  given  by  the 
eigenfunction  |R>  of  the  nuclear  position  operator.  To  this  end 
we  use  the  resolutions  of  the  identity  operator  in  nuclear  space 


S dR  |  RXR  |  =  1 

and  in  the  electronic  space 
/dij)  1  |  =  1  . 

Here 


X dd>  a  (N+l  )  ! 


N  d(ReXr,)d(Iaixn)  ) 


n=l 


) 


4 . 3 


4 . 4 


(4.5) 


where  the  prime  denotes  integration  under  the  restriction 
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that 


The  identity  ( 4 . 4 ) - { 4 . 5 )  may  be  easily  proven  by  showing 
<r.  |  |  i(i><$  |  m>  =  (  N+l )  !  S  1  (  dx-^dy./n)  ( :<n~  iyn  )  (  x^eiy^ )  =  5^. 


The  path  integral  can  be  constructed  in  the  usual  way 


K  {  R’4i  ,  R  4)  |  t ) 

o  '  o 


<RL*L 


j  =  l 


-  f  Kit 

e  71  I  R  4i  > 
o  o 


f  L~1 

FT  jdR  /dm. 

n=i  j  j 


n 


j=i 


,  .  -lHA t/n  .  _ 

<  R  .  '!>  .  i  e  R  .  ,0). 

J  J  J~1  J- 


4.6) 


where  At=t/L,  4.  =  41  ( t  .  )  ,  t.=t  +  j  A  t  and  where  R  and  4>r  stand  for 

R  and  >4/  at  time  t.  Focussing  now  on  the  matrix  element 

<R  .  4)  .  I  eHA  t /h  |  R  .  ,  41  •  ,  >  we  evaluate  its  electronic  cart  first, 
j  j  J  ~  1 - 1 

keeping  R  and  ?  in  the  Hamiltonian  (4.1)  as  parameters.  For 
completeness  we  outline  the  procedure  in  Appendix  A.  The  result 
is 


—  i  H  A  t  /  n 

<R  .ill  .  i  e  '  R  .  '!>  .  ,  >  = 

J  J  ’  3-1 ■ J-l 


4.?) 


.-‘VV«,i.‘  6%  -  V*”  V 

J  J 


» 
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Here  P  and  R  are  operators  in  the  space  spanned  by  the  nuclear 
wave  functions,  and  <  ^  >  is  given  in  terms  of  the  coefficients 

x  by 


=  <(Z  xn(  tj  )  |n>  )  |  (~  Z  xn(tj)|n)> 


(4.8) 


"  2  Xn  (t;j)xn(tj]  . 


V.(R)  is 


V  .  (R)  =  Z|xn(t  )  |2H  (R)+Z  Z  x*  (t  )xm(t.)H  (R) 

1  n  n  j  nn  n  m  n  jmjnm 


(4.9) 


The  nuclear  matrix  element  in  (4.7)  is  now  seen  to  take  a 
form  which  is  normally  obtained  for  a  quantum  particle  moving 
under  the  influence  of  a  potential  V(R).  Evaluating  the  R  matrix 
element  in  the  standard  way  and  inserting  it  into  (4.6),  then 
taking  the  L  -*>  °°  limit,  yields 


ROi’ 


K {  R(*i ,  Rq4»0  !  t )  =J  DiiiDReiS(R'J)  ^o^o1  t 


4.10' 


R  41 
oTo 


with 


S(R<|i,Romo|  t)  =  J  dt 1  L  (  R  ( t  ’  )  ,  R  ( t 1  )  ,x(t)  ,  x  (  t  '  )  ,t'  )  (4.11) 


and  the  "Lagrangian"  L  given  by 

,  2. 


-  *  Z 
n  n 


■n 


-(H) 

nn 


Z  Z  x 
n  m 


n  m  nm 


:r>>  : 


(4.12' 


the  symbols  Dv  and  DR  are 
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L-I  N  d(ReXl(tJ) )d(Imx1{tj) 

□  4>  =  lim  1  [  {  ( N+l )  '1  [  T - 

L-~»  ;j  =  l  1  =  1 


4.13; 


L  L-l 


DR  H  11 


im  {(J 


2TTih(  t/L; 


dR  ( t  )  } 

j 


(4.14) 


Equation  (4.10)  Is  an  exact  representation  of  the  propagator  in 

terns  of  a  path  integral  over  all  the  paths  which  lead  from  the 

state  R  ,  t|i  =  Ex  (  *  )  |n>  at  time  t  to  the  final  state  R  and  iji  = 
o  o  n  o  o 

2xn(t)|n>,  at  time  t.  Having  found  this  form  it  is  of  interest 

to  look  for  the  optimal  path,  the  analog  of  the  classical  limit, 

as  an  approximation  to  the  path  integral.  For  this  purpose 

define  X  =  Rex  ,  Y  =  Imx  and  replace  in  Equation  (4.12)  x  = 
n  n  n  n  ^  ^  *n 

* 

X  +  iY  and  x  =  X  -  iY  everywhere.  To  satisfy  the 
n  n  ^n  n  n  1  1 

2  2  2 
restriction  I]x  I  =1  we  may  add  to  L  a  term  XI  ( X  +Y  )  where  X  is 
n  n  n 

a  lagrange  multiplier,  but  it  turns  out  that  this  does  not  change 
the  result.  The  optimal  path  is  the  solution  of  the  Euier- 
Lagrange  equations  of  motion  generated  by  minimizing  the  action 
given  by  Eqs .  (4.11-4.12): 


dt  3 R 


dt  3X 


3  L  d  3  L  „  ,  ,  , 

—  -  at  vr  =  0  (al1 

n  n 


(4.15) 


These  yield  the  following  Equations  of  motion 


.  3  V  (  R  ,  {  x  _  '  ) 
M  3  R 


(4.16 
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x„  = 


K[Hnn< 


H 


nl< 


(R)x,] 


;  n=  i ; 


with  the  effective  potential 


V(R, {xn> ) 


=  Z 
n 


H  <R)  !  x 


nn 


n 


Z  Z 
n  m 


H  (  R  )  x  X 
nm'  ^nAm 


(4.18) 


To  obtain  Equation  (4.17)  we  have  recombined  the  equations 
obtained  for  and  Y  .  Equations  ( 4 . 1 6 ) - ( 4 . 1 8 )  are  identical  to 
the  equations  of  motion  derived  in  Section  III. 


The  following  comments  should  be  made  with  respect  to  this 

derivation:  (a)  The  same  results  are  obtained  If  instead  of 

* 

replacing  xn  by  xn+— n  we  regard  xn  and  xn  as  independent 
variables  and  use  the  Euler-Lagrange  equations  in  the  form  3L/3xn 
-  d/dt ( 3 L/ 3 xn ) =0 . 

(b)  In  the  absence  of  coupling  the  nuclear  motion  (i.e.  if 

we  were  evaluating  the  time  evolution  of  an  N-level  system)  the 

resulting  EO.t  (4.17)  is  just  the  Schrodinger  equation  for  this 

system.  Obviously  the  optimal  path  is  over  determined  by  the 

requirements  that  (  t Q  =  !-P 0  and  0  (  t )  =4i  are  given,  and  generally  we 

will  not  find  a  path  that  will  satisfy  both  these  requirements 

and  Eq  (4.17)  .  The  result  may  be  still  interpreted  in  the 

following  way.  The  probability  amplitude  to  be  in  a  final  state 

starting  from  an  initial  state  ii  may  be  written  in  the  form 

o 

do  <w  _  ,  li  >X  (  ii  ,  lu  I  t )  .  The  optimization  procedure  replaces  this 
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integral  bv  the  maximum  value  of  the  integrand  which  is  obtained 


it  >D=I\  |  n>  is  the  solution  of  Eq.(4.17). 


The  procedure  described  above  may  be  used  to  obtain 
equations  cf  motion  for  more  general  si  tua  t  ions  i  nvo  1  vi  r.g  coupled 
quantum- classical  systems.  For  example,  the  Hamiltonian  (4.1) 
can  be  supplemented  with  coupling  terms  arising  from  the  effect 
of  the  nuclear  kinetic  energy  operator  on  the  electronic  states, 
leading  to 


H  (R)  | n>  <m |  + 
nm  M 


W  ( R ) | n> <m | 
nm 


n ,  m=  1 
_  2 


n ,  ra=  1 


(4.  IS) 


+  ~  I | n>  <n | 
n 


w  1 1  n 


2  2 

H  (  R )  =  H)  (R)6  mV  (R)-t-  i'd3rp*  (  r  ,  R  )  - — r  o  (r.R) 

nm  n  nm  nm  2.  m  n  „  —  *  m 

o  R 


(4.20) 


anc 


W  ( R )  =  - i h  S d “ r  6  (r.R)  tt  b  ( r , R ) 

nm  n  «R  m 


(4.21' 


Here  b„(r,R)  a  | n>  is  used  to  denote  the  explicit  dependence  of 
the  electronic  wave  function  on  the  electronic  (r)  and  nuclear 
(R)  coordinates.  The  procedure  described  above  may  be  used  to 
obtain  the  propagator  in  the  form 


^  i 


J  .22) 
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;econd  mechanism  consists  in  direct  energy  transfer  from  the 
noving  ion  to  the  electrons  in  the  metal,  and  it  is  discussed 
nere  in  detail.  This  also  has  two  distinct  aspects.  One 
corresponds  to  electron-hole  pair  excitations  caused  by  the 
sudden  creation  of  the  ion  through  electron  transfer  and  it  was 

discussed  previously  in  a  different  physical  context  by  Gadzuk 

4  3  ••  4Q 

and  Metiu  and  Schonhammer  and  Gunnarson.  The  other 

corresponds  to  electron-hole  pair  excitations  by  the  ionic  motion 

3  7  3  4 

and  was  considered,  for  example,  by  Suhl  etl  a  1  .  '  **  and 

Korzeniewski  et.al.04 

In  the  present  section  we  consider  the  effect  of  the  energy 
transfer,  to  excite  electron-hole  pairs,  on  the  motion  of  the 
ion.  In  the  spirit  of  the  Langevin  approach  to  many  body 
dynamics  we  woulo  like  to  produce  an  equation  which  does  not  make 
explicit  reference  to  the  detailed  electron-hole  pair  dynamics. 
Since  the  motion  of  the  ion  drives  the  excitation  of  the 
electrons,  it  must  provide  the  excitation  energy:  therefore,  the 
desired  equation  must  contain  a  friction  term.  Furthermore,  the 
interaction  of  the  ion  with  the  polarization  fields  caused  by  the 
thermal  fluctuations  of  the  electron  in  the  metal  must  be  taken 
into  account  through  the  presence  of  a  random  force  in  the 
equation  of  motion  of  the  ion.  Since  the  metal  screens  the  ion 
field  we  can  assume  that  the  screened  interaction  causes  a  small 
perturbation  in  the  motion  of  the  electrons  in  the  metal 
Therefore  we  expect  t.ne  random  force  ana  the  friction  kernel. 


not  only  the  nuclear 
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The  coupling  to  the  lattice  affects 

motion  but  also  the  transition  amplitudes  c(t)  since  H  ,(R,Y,X  ) 

i  i  j  o 

are  now  function  of  the  stochastic  variable  Y.  Both  the  crossing 
point  and  the  coupling  strength  are  random  variables.  The 
observable  quantities  of  the  theory  must  be  computed  by  running 
many  stochastic  trajectories  and  averaging  the  quantity  of 
interest  over  them. 

VI. 2.  Energy  Loss  Due  to  Electron-Hole  Pair  Excitations. 

Recently  there  has  been  a  lively  interest  in  the  manner  in 
which  the  excitation  of  electron-hole  pairs  in  the  surface 
influences  the  dynamics  of  adsorbates . 3 ■ 4  *  9 • 10  ■  12  ■  21 ■ 35  55  Two 
types  of  mechanisms  are  possible  in  the  problem  of  interest  here. 
In  the  first,  the  electron  jumping  from  the  atom  into  the  metal 
interacts  with  the  electrons  already  there  and  it  is 
inelastically  scattered.  If  the  interaction  partner  is  below  the 
Fermi  level  the  interaction  can  excite  an  electron-hole  pair. 

This  is  similar  to  a  shake-up  process  and  it  could  be  represented 
by  giving  each  ionic  state  a  width  which  corresponds  to  the  rate 
of  energy  loss  from  the  transferred  electron  to  the  electrons 
below  the  Fermi  level.  The  transferred  electron  can  also 
interact  with  the  electrons  thermally  excited  above  the  Fermi 
level  and  this  leads  to  either  energy  loss  or  gain.  This  can 
also  be  represented  by  giving  the  ionic  state  a  width.  The 


•>  2 
Md  Y/ dt' 


=  -  X,  3  H  (  R  ,  Y  .  X  )  /  3  Y 
1  11  o 


2Re  x.X,  3  H  (R,Y . X  ) / 3  Y 
12  12  o 


(6.3) 


•  i  xx 

M-l  *  " 


Here  X  ^  is  the  force  constant  coupling  X(i  to  Y.  The  equation  of 

motion  for  the  secondary  atoms  is  that  of  the  harmonic  lattice  in 

the  absence  of  the  ion.  X  represents  the  coordinates  of  the 

o 

secondary  atoms  in  their  equilbrium  postions. 


Following  the  Adelman-Doll  procedure  we  can  eliminate  the 
secondary  lattice  atoms  to  obtain  a  Langevin  equation  for  the 
primary  atom : 


MdV'dt2  =  -  Xj  ;  2(  3H11/3Y)-2Re  x*  x2(3HJ2  3Y) 


/  y( t-t ’  ) Y( t  '  )dt 1  -  F ( t ) 


The  friction  kernel  is  proportional  to  the  linear  response 
function  of  the  lattice  in  the  absence  of  the  primary  atom.  The 
random  force  F(t)  is  Gaussian  and  its  correlation  function  is  the 
Green's  function  of  the  secondary  lattice  in  the  absence  of  the 
primary  atom.  The  Green's  function  and  the  linear  response 
function  are  related  by  the  fluctuation-dissipation  theorem.  The 
Langevin  equation  can  be  solved  efficiently  by  modeling  the 
friction  and  generating  Fit)  on  the  computer  as  shown  by  Shugard 


e  t  .  a  1 


under  the  influence  of  the  effective  potential 


2  2  , 

V ( R , X ; t )  -  Z  Z  X,(t)x,(t)  H(R.X)  .  (6.1) 

i = 1  j=l  J  J 

We  assume  here  that  the  neutral  curve,  which  we  denote  H 22 ,  is 
independent  of  X.  In  other  words,  the  energy  lost  by  the  neutral 
to  phonons  is  disregarded  since  in  this  context  it  does  not  have 
marked  effects  on  sticking.  It  can  however  be  easily  included, 
if  neces  sary . 

In  order  to  derive  a  Langevin  equation  for  the  present 

3  3 

situation  we  follow  Adelman  and  Doll  and  divide  the  lattice 
atoms  into  a  primary  zone  which  suffers  the  brunt  of  the 
collision  with  ion,  and  a  secondary  zone  which  interacts  with  the 
primary  atoms  only.  For  simplicity  we  take  only  one  primary  atom 
and  denote  its  coordinate  by  Y(t);  the  other  lattice  atom  form 
the  secondary  zone  and  their  coordinates  are  denoted  either  by 
X  .  u=l....N,when  we  need  to  specify  all  of  them,  or  by  X.  when 
they  are  denoted  collectively.  The  equation  of  motion  for  the 
ion  is 

2 

>2  * 

md“R/dt  =-I  xx.  3H..(R(t),Y(t).X)/3R,  (6.2) 

.  ,  l  J  l  J  o 

l  ,  j  =  1 

The  equation  of  motion  for  the  primary  zone  atom  is 


VI 


THE  INCLUSION  OF  PHONON  AND  ELECTRON-HOLE 
EXCITATION  AS  ENERGY  LOSS  CHANNELS. 


As  we  have  discussed  in  Section  II  if  the  kinetic  energy  K 
of  the  incident  neutral  is  higher  than  the  asymptotic  energy 
mismatch  A  V  ( » )  (i.e.  A  =K-AV  ( °°  )  >  0  )  the  ion  formed  by  charge 

CX  Ct  CX. 

transfer  sticks  to  the  surface  only  if  it  loses  (to  phonons  or 

electron  hole  pairs)  an  amount  of  energy  larger  than  A  In  what 

a 

follows  we  outline  a  curve  crossing  theory  which  include  such 
energy  loss  processes.  Both  consist  of  deriving  Langevin 
equations  in  which  the-  action  of  phonons  or  electron-hole  pairs 
generates  friction  and  random  forces  in  the  mean  trajectory 
equation. 

VI. l  The  Inclusion  of  Phonons. 

We  consider  here  the  two  state  problem  discussed  in  Section 
III. 3  and  include  in  both  the  wave  function  and  the  electronic 
Hamiltonian  the  coordinates  of  all  the  lattice  atoms.  For 
simplicity  we  denote  all  these  coordinates  by  X.  and  consider  the 
one  dimensional  case  only.  The  extension  to  three  dimensions  and 
more  than  two  electronic  states  is  straightforward. 

By  repeating  the  derivations  presented  in  Section  III. 3  we 
obtain  equations  identical  to  (3.35)  and  (3.36)  in  which  the 
matrix  elements  H  .  ^  depend  on  R(t)  and  X!t).  We  must  now  specify 
a  procedure  that  gives  the  equation  of  motion  for  R ( t )  and  X ( t  ) 
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and  where  the  prime(')  denotes  derivative  with  respect  to  R.  In 
1'^  this  derivative  is  obtained  by  replacing  ( R ) ^ ( R )  by  its 
R  derivative  with  H^(R)  kept  fixed).  The  set  of  equations 
( 5 . 22 ) - ( 5 . 23 )  may  be  now  integrated  numerically  if  a  model  for 
the  coupling  matrix  elements  is  constructed. 
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C1  = 


-  ^  /  dT£21 ( T  * 


(5.22a) 


i[K21(R)-Z21(H)]e 


X  dT£12(T) 

Cl+  i  S22(R)C2 


( 5 . 22b) 


raR  =  -{  (H^(R)-tr  q(R)  )  (cj2  +  (  H^2  (  R) -R-  q  (  R)  )  |  c2  |  2 


-  H^q(R)  +  2Re [Z2 ( R) \c^ 


(5.23) 


,  h  l  dT£2i(T)  , 
-(H^_2(R)-Z^2(R)-2^1(R)  )e  C1C2]} 


where 


£  2 1  =  H2 2  “  H11 


(5.24) 


Zkl  Dkl  irki 


Oc,l  =  1,2) 


(5.25) 


H.  g. ( R)  H  ,(R) 

Dkl(R)  =  ppfdE°(E}  hqo(R)+e-h11(R) 


5.26) 


R)  =  TTXdEp(E)HlcE(R)HEI(R)6(Hoo(R)+E-H11(R)  ) 


P 
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to  those  described  above  where  we  now  have  two  discrete  states  i 
(neutral)  and  2  (surface  state)  coupled  to  the  continuum  of  ionic 
states.  The  neutral  state  coupling  to  the  surface  state  leads  to 
the  ionization  of  the  neutral.  The  electron  transferred  into  the 
surface  state  can  be  recaptured  by  the  ion  (the  coupling  is  H2g) 
leading  to  neutralization.  The  mean  potential  for  this  case  is 


V(R(t)  ) 


2  2  2 

=  I  Z  H  ,(R(t)  )x?x,+  Z  J,dep(e)[H.{R(t)  )x*x 
i  =  l  j  =  l  J  j  =  l  J  " 

(5.19) 

+  Hej{R(t)  )x*xj]+;dep(e)Hee(R(t)  )  |  x£  I  2 


The  classical  equation  for  R(t)  is,  as  before, 


MR  =  -3V ( R ) / 3R  (5.20) 

while  the  equations  for  the  amplitudes  x  are 

X.“  -  ^(H:i(R(t) )x:+H12(R(t) )x2+/dep(e)Hle(R(t) )x£) 

i2  «  -  S  (H22(R(t)  )x2+H21(R(t)  )x1^J'deo(e)H2e(R(t)  )x£) 

Xe  -  -  £  (HeeR(t))xe+H6l(R(t) )Xl+He2(R(t) )x2) 

(5.21) 

Proceedings  along  lines  similar  to  those  described  below  Ec. 
(5.4)  leads  to  the  set  of  equations 
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+®  dE  p(E) | H 
)  =  P?  S 


E+H  ( t ) -H . 
-oo  oo  '  11 


11 


and 


(t)  =  tt  ;  dE  p(E)  |H1Er6(E+Hoo(t)-H,,(t 


(5.12; 


If  the  continuum  states  form  a  narrow  band  (which  means  that  o{s) 
is  aero  outside  a  narrow  energy  range)  these  equations  will  be 
modified  by  "edge  effects".  A  very  narrow  ionic  band  influences 
the  neutral  state  just  like  a  single  ionic  state. 


The  - tlculations  car 

ried  out 

above 

e  f  f  ec 

tive  potential  given 

by  Eq.  ( 

5.4)  i 

and  r 

( t )  alone .  The  last 

term  in 

Eq .  (5 

permit  us 
terms  of 
4 )  can  be 


to  rewrite  the 
'c.ft);2,  A(t) 
written  as 


W e  obta; n 

this  by  observing 

that 

~EE  1  "  '  r'oo  '  )  ~~ 

and 

that  we 

r.eed  only 

the  gradient  cf  H_ 

,  (  to 

compute  the  force 

)  so 

we  can 

t  a  .<  e  H _ ;  t 

;  -  J>  (  t  '  !  .  A 1 

so  we 

use  the  normal iza 

t  ion 

condi ti 

•  c .  (  t ) 


(5.14) 


The  middle  two  terms  in  Eq.  (5,4)  can  be  rewritten  by 
expressing  xc  in  terms  of  x-  and  using  the  approximations 
presented  under  the  equation  (5.7);  the  result  is  that  the 

n 

terms  are  equal  to  - 2*. ( t )  f c .  •  t }  '  "  . 


wo 
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is  a  time  independent  energy  giving  the  difference  between  the 

parallel  curves  H  and  H  .  If  we  assume  that  the  "amnlitude" 

LL  OO 

H, _ ( t ' ) c  ( t ' )  varies  with  t'  on  a  time  scale  much  slower  that  the 
It  1 

phase  (i(t)-n(t')  +  E(t-t'),  than  we  can  replace  H  (t')c  ( t '  ) 
with  H  w(t)c  (t)  and  remove  it  from  the  integral  over  t'.  After 
that  we  can  perform  the  integral  over  t'  by  inserting  a  factor 
exp(+r|t)  ti>0,  to  insure  that  the  integral  converges.  We  get 


t  t 1 

J  dt'  exp { - ( i /h )  !  [ f ( x ) -e+ig ]dt } 
—  00  —00 


(5.3) 


=  -  (fr/i  )  (  f  ( t ) -e  +  iTi )  exp  {  -  (  i  /h )  f  [  f  (  t  ) -e  +  irj  ]  dx  ; 

—  00 


in  deriving  this  we  assume  that  f(r)  (which  in  our  case  is 
H1 i ( r ) -Hq0 ( T ) +E )  can  be  replaced  by  f(t)  and  be  treated  as  a 
nearly  constant  quantity.  After  performing  these  calculations  we 
obtain : 


ihc  x  (  t )  =  [A  (  t )  -il"(  t )  ]  c  :  ( t ) 


(5.9) 


with 

A(t) 

r(  t) 


Re  > 

■  =  lim 

dSp (E) 

H..P(t) 

X  w 

/  rr*  0 

Im  . 

J 

CE+Hoo(t)-H11(t)  +  ir1] 


(5.10) 


If  the  integration  interval  is  practically  infinite  the  last 
equation  leads  to 


disregarded.  The  effective  classical  interaction  is  given  by 


V(R(t))  =  H11{R(t))Ix1(t>!2  +  Jde  p(e) [Hlf(R(t) )x1(t)xe(t) 

+  H el(S(t))x  (t)x.(t)]  +  Sde  p(e)H  (R(t))|Xp(t)|2  . 

(5.4) 

The  "classical"  trajectory  can  be  computed  by  using  Hamilton 

2 

equations  with  the  Hamiltonian  P(t)  /2m+V(R(t)). 

We  can  simplify  these  equations  by  following  a  standard 
procedure  used  in  situations  when  a  discrete  state  interacts  with 
a  continuum.  We  start  by  integrating  exactly  Eq .  (5.3)  to 

express  \£  as  a  function  x1:  we  introduce  this  equation  for  x£  in 
Eq .  (5.2)  and  obtain  a  closed  equation  for  x1 

t 

ihc^t)  =  J*  dt'  G  ( t  •,  t '  )  c  ( t '  )  (5.5) 

—oo 

with 

G ( t ,  t '  )  s  ( -i/fr)  /dEp(E)H  (t)H, „ ( t  '  ) 
exp{  - ( i/h) [p( t ) -u( t ' ) ] )  exp { - ( i/h) E( t-t '  ) )  , 


and 


h(  t ) 
x  1  ( t ) 


t 


s  ;  [Hll(T)-Hoo(T)]dT 
-00 

V. 

=  c1(t)  exp{ - ( i/h)  X 


H11(T)dT) 


(5.6) 

(5.7) 


Here  we  use  the  fact  that  the  ionic  curves  are  parallel  to  each 

other.  We  pick  one  of  them  as  a  reference  state  and  denote  it  by 

H  (R(t))  sh  (t),  and  write  the  others  as  H  (R(t))+E,  where  E 
oo  oo  oo 
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V.  CROSSING  OF  A  DISCRETE  STATE  BY  A  CONTINUUM 


As  we  discussed  in  Section  II  if  a  neutral  atom  approaches 
the  surface  it  can  be  ionized  by  crossing  to  an  infinite  number 
of  states  differing  from  each  other  through  the  final  state  of 
the  transferred  electron,  the  position  of  the  crossing  point  and 
the  amount  of  asymptotic  energy  mismatch  tc  be  compensated.  To 
treat  this  case  of  a  neutral  curve  crossing  a  continuum  of  ionic 
curves,  by  the  methodology  used  in  Section  III.  3  for  cne  ionic 
state,  we  use  for  the  amplitude  of  the  total  wave  function  Ec . 
(3.10)  the  form 


A  (  X  ,  R ;  t )  =  x1(R)4>1(x,R)  +  Jde  p ( e ) x£ ( H ; t ) 0£ ( x , R )  (5.1) 


The  procedure  of  Section  III. 3  applied  to  the  wave  function 
(3.10),  with  the  amplitude  (5.1),  leads  to 


ihx,  (  t )  =  H.  (?.(  t)  )x.  (  t)  +  /as  p(s)H._  (R(t))x.(t) 

I  11  x  1  £  t 


•*  \  /  R 


inxe(t)  =  H££  (R(t))x£!t)  -  H£l(R(t))x:  ■ 


v  -  ■  J  / 


We  neglect  the  matrix  elements  H  ,  ( e^e ' ) ,  since  the  continuum 


states  do  r.ct  interact  with  each  other.  In  the  present  context 
this  means  that  the  motion  of  the  ion  does  net  induce  transitions 
of  the  transferred  electron  from  e  to  s' ,  and  also  that  electron 


r*.  c  is  Zj  L  z"  0*<cit3ti.on  of  rn  0  tr  s  2  fe  v  th°  i  c  n  1  c  xotiicr*  is 


I 

I 
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with 


S  =  fdt '  [i<ip(it>> 
o 


2 

K  lP*~  Im  "  fi  w(r,40-u(r,<m  )  ] 


14.23) 


where 


W(R,4>) 

=  Z 

Z 

X*X  w  (R) 
n  m  nm 

n 

m 

*  — 

U(R,i|0 

=  Z 

Z 

X  X  H  (  R  ) 

n  m  nm 

n  m 


(4.24) 


Here  we  have  used  the  phase  space  (R,P)  representation  of  the 
nuclear  part  of  the  path  integral.  The  ROM's  for  the  optimal 
path,  obtained  from  the  extremum  of  the  action  S  given  by  Eq . 
( 22 ) ,  are 


•  _  _  P  3W(R,H) )  _  3  U  (  R  ,  ij) ) 

?  M  3R  3R 

R  =  ^  (  P  +  W  (  R ,  i|i  )  ) 


V 

n 


r>  rn 

•  iitt 


( R ) 


o 

^  w 

M  nm 


(  R )  ]  x 

m 


The  implications  of  these  equations- will  be  explored  elsewhere. 


appearing  in  the  equation  of  motion  of  the  ion,  to  be  related 


▼ 
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f  e  (  r  .  r  ;  u  )  £  (r')dr’  =  o 

q  q 


(6.7) 


are  used  to  obtain  V  (R(t))  from  the  equation 


V  ( R  (  t )  )  -  e2Z  J*  |  "r  -  R(t)1  1g  (r)dr  .  (6.8) 

H  H 

V  (R)  is  the  screened  ion-electron  coupling.  Here  e  and  Ze  are 

q 

electron  and  ion  charges. 


We  now  wish  to  obtain  the  equations  of  motion  in  which  the 
coordinate  R(t),  its  conjugate  momentum  P(t)  and  q  and  r\~  are 

q  q 

classical  variables  moving  in  a  mean  potential  analogous  to  the 
one  obtained  in  the  previous  sections.  The  procedure  used  in  the 
previous  sections  leads  .to  the  following  classical  Hamiltonian 


K^R.P.q.q  -.XjXg) 


p  JR 

—  *  Ihwriri 
2ni  q  q  q 

q  s 


x,:2Hnm-  x2:2H22(R>-x;x2H12m-x,X2H2i(RI 

-  x, : 2  i  <VR,VVR)V 


(6.9) 


o : 


* 

Note  Chat  q  .  the  complex  conjugate  of  appears  in  place 
the  operator  n ‘  .  in  accord  with  the  classical  equivalent  of 


the  boson  creation  operator 


From  this  Hamiltonian  the  equations  of  motion  in  the  mean 
trajectory  approximation  are  obtained  using  Hamilton's  equations 


» 


R 


3  H  /  3  P 
c 


P/  m 


(6.10) 


P  =  -3H  /  3R  =  -  lx.  !  2[3H  /3R+I{  (  3  V  /3R)rj*-*-(  3V*/3R)ri  } 

C  1  11Q  Q  Q  Q 


3H22/3R- 

2Re(XjX23H12/3R) 

(6.11) 

3  H  /3ti* 
c  'q 

■■‘vr  s 

(6.12) 

3  H  /3t| 
c  'q 

*  K  VR"*ii2 

(6.13) 

In  addition,  the  quantum  equations  of  motion  for  the  amplitudes 
Xj  and  x2  are  obtained  in  the  form 


i‘hx1  =  [H11(R)+E{Vq(R)T1ci+Vci(R)r1c?)  ]Xj  -  H12fR)x2  (6.14) 

i*V*2  =  H22(R)X2  *  H21(R)Xj  (6.15) 

The  equations  (6.12)  and  (6.13)  can  be  formally  integrated 
to  give 


t  ,  iw  (  t  '  -  t  ) 

T)q  (  t )  =  Tiq(t)  -  X  (i/h)  Vq  (  R  (  t  '  )  )  :  x  x  (  t  '  )  e  q  dt' 

(6.16) 

o  o  -  i  u)  t 

where  ti  (t)  =  t|  (o)e  q  is  the  value  of  t)  (  t  )  in  the  absence  of 

q  q  q 

the  ion.  The  initial  value  ri°(o)  may  be  regarded  as  a  stochastic 
variable  given  by  the  (classical)  thermal  distribution 
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TT2k2T 


exp { -  (  4hu^/kT)  !  j  2  } 


(6.17) 


Note  that  Planck's  constant  appears  in  this  classical 

distribution  because  it  is  used  to  construct  the  dimensionless 

* 

quantities  r>^  and  from  the  momentum  and  the  coordinate  of  the 
oscillator . 


We  summarize  these  results  by  outlining  how  they  should  be 

used  in  a  calculation  in  which  all  practical  concerns  (i.e. 

computer  cost)  are  disregarded.  Let  us  assume  that  we-  want  to 

calculate  the  thermal  average  of  a  quantity  A(t)  = 

A( ( t ) , x2 ( t ) , R( t ) ) .  This  is  the  expectation  value  of  the 

dynamic  quantity  A  (i.e.  A  ( x,  ( t )  .  x2  ( t )  ,  R(t))  2  <H>  (  t )  |  A  |  41  ( t )  >  ) 

and  can,  in  gene'al,  depend  on  x 1 , X 2  and  3.  To  start  the 

calculation  we  generate  by  Monte-Carlo  the  values  of  r\  (o)  and 

cr 

* 

n  (o)  for  all  the  bosons  (i.e.  q=l,2,...N)  reauired  for  a 
c 

realistic  description  of  the  electron  excitations  of  the  gas 

* 

svstem.  We  denote  the  set  {*n  (o),r|  ( o )  > _ .  v  by  p.  We  can 

q  q  q- * 

now  solve  simultaneously  the  Eq  (6.16b)  for  -i  (  t }  and  n  ( t ) ,  the 

Cl  c 

Ecs .  (6.14)  and  (6.15)  for  x,(t)  and  x-(f),  and  the  equations 

*  ±  * 


(6.10)  and  (6.11)  for  R( 


This  orovides  the  cuantities 


?.(t),x,  (t)  and  x^[t)  ,  for  any  desired  time.  We  can  therefore 
compute  A  (t),  corresponding  to  the  set  p  initial  condition. 
Using  important  sampling  Monte-Carlo  we  can  generate  mere  ir.it: 
conditions  and  repeat  the  calculation  outlined  above.  The 
thermal  average  <A>  of  A  is  given  by 
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<A>  =  Z  A  (X:(t) ,x2(t) ,R(t) )  . 

P  P 


6.18) 


While  this  calculation  illustrates  the  structure  of  the 
theory,  it  is  rather  expensive.  A  simpler  procedure  is  to  write 
(by  using  Eq.  (6.16)) 


S(v  (R,^*(t)  +  V(R)V)  =  ?(t)  +  /  r(t-x)  |x,(t)  |  2dx 
q  q  q  q  q  o 


with 


(6.19) 


S(t)  =  2Re  Z  V  (t)r,®(t) 

q  q 


r(t-T)  =  Re ( 2/ih)  Z  V  (t)V  (T)exp(iw  (t-T)] 

q  q  q  q 


6.20) 


We  can  now  treat  f(t)  as  a  Gaussian  variable  at  each  time  t  with 
the  correlation  function 


<*(t)s(T)>  =  Re  Z  V  ( t ) V  (t) (kT/hu  )exp(iu  (t-T) ] 

q  q  q  q  q  (6.21 


We  can  model  the  sum  above  by  using  a  simple,  but  appropriate, 
function  of  t  and  t;  this  can  save  substantial  computer  time, 
since  we  need  to  generate  one  Gaussian  variable,  i.e.  f(t),  for 
each  time  point.  In  making  such  models  we  must  assure  that  the 
exact  relationship  between  y(t-x)  and  <?(t)§(T)>  is  preserved; 


otherwise  fluctuation-dissipation  theorem  is  violated. 
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VII  DISCUSSION 

The  main  thrust  of  this  paper  is  that, the  theory  of  low 

energy  collision  phenomena  which  can  be  described  as  two  or  many 

state  problems,  may  have  to  discard  the  trajectory  approximation. 

The  reason  for  this  is  rather  general  and  the  problem  is  present 

whenever  we  deal  with  two  sets  of  degrees  of  freedom  x  and  R  and 

we  must  treat  x  quantum  mechanically,  but  intend  to  treat  R 
29 

classically.  If  we  replace  R  in  the  Hamiltonian  with  a 
trajectory  R(t),  obtained  by  following  some  classical 
prescription,  we  generate  a  time  dependent  Hamiltonian  for  the 
quantum  variables  x;  this  forces  the  quantum  subsystem  to  undergo 
transitions  fron  the  initial  state  ji>  to  some  set  of  final 
states  | f > ,  f=l , 2  ...  The  system  conserves  energy  only  if  the 
energy  loss  (or  gain)  along  the  trajectory  R(t)  equals  E*-Z,. 
Since  the  latter  is  different  for  each  final  state,  energy 
conservation  requires  a  different  trajectory  R.  e ( t )  for  each 
transition  i  -*  f .  Therefore  the  force  ?.  *(R{  e  ( t )  )  appearing  in 
the"Newton"  equation  for  R.  *(t), 
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simplicity  we  assume  that  the  interaction  between  R  and  x  is  non¬ 
zero  only  during  the  collision;  this  is  true  for  the  charge  and 
energy  transfer  problems  considered  here) .  Since  energy 
conservation  fixes  the  final  state  of  the  classical  degrees  of 
freedom  we  expect  F.  ^(R.  -  ( t ) )  to  be  a  functional  of  the 

trajectory  and  the  Eq.  (7.1)  to.be  an  integro-dif f erential 
equation.29,55'5®  Furthermore,  due  to  the  nature  of  quantum 
mechanics  each  transition  i  -*•  f  for  the  quantum  subsystem  has  a 
probability  W^_^(R^  ^(t))  which  is  a  functional  of  the 
appropriate  trajectory  R^  f(t).  Thus  a  complete  description  of 
all  the  events  possible  in  such  a  system  consists  of  a  catalog  of 
all  the  transitions  i  -*  f  that  do  not  violate  conservation  laws, 
a  set  of  integro-dif f erential  "Newton's"  equations  of  the  form 
(7.1)  which  give  a  trajectory  for  each  transition,  and  a  set 
consisting  of  one  transition  probability  per  transition. 

These  general  considerations  do  not,  however,  tell  us  how 
to  derive  these  equations.  The  obvious  strategy  is  to  formulate 
the  problem  fully  quantum  mechanically  and  then  to  take  the 
classical  limit  for  R.  We  have  explored  three  such  methods:  one 
using  the  stationary  phase  approximation  in  the  path  integral 
expression  for  the  transition  amplitude55'57:  a  second,  using  of 
wave-packets  to  generate  classical  like  equations  for  the 

3 1 

position  of  the  center  of  the  packet  and  its  group  velocity:” 
and  a  third,  which  uses  the  eikonal  method.  The  classical 
approximation  for  R(t)  can  be  carried  out  within  each  method  at 
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various  levels,  generating  different  "classical"  theories  for 
R(t).  The  mean  trajectory  approximation  ( MTA )  presented  here  is 
the  simplest  of  them. 

Within  a  path  integral  representation  the  MTA  is  obtained 
by  using  a  coherent  state  representation  for  the  electronic 
degrees  of  freedom  and  by  obtaining  equations  of  motion  for  all 

C.  £? 

the  variables  by  using  the  stationary  phase  approximation. 

Thus  the  electronic  amplitudes  are  treated  as  classical  fields. 

In  spite  of  this  we  get  for  them  a  time  dependent  Schrodinger 

equation  driven  by  an  effective,  .mean  potential.  If  we  were  to 

take  a  classical  limit  far  the  nuclear  coordinates  only,  then  we 

5  7 

would  have  obtained  one  integro-dif f erential  equation  for  the 
trajectory  R^(t),  depending  on  the  initial  and  final  electronic 
state . 

Within  the  wave  packet  formalism  MTA  is  obtained  by  faking 

one  Gaussian  packet  for  both  electronic  states.  A  refined 
5  8 

version  takes  t-wo  packets,  one  for  each  electronic  state,  and 
generates  two  "classical"  equations,  for  two  trajectories  which 
jointly  drive  the  electronic  degrees  of  freedom. 

Finally  the  eikonal  method  used  here  can  be  refined  by 

5  8 

defining  two  eikonals,  one  for  each  electronic  state. 

This  brief  survey  of  various  possibilities  serves  to 


pinpoint  the  limitations  of  the' mean  trajectory  approximation. 
Since  it  generates  only  one  trajectory  it  satisfies  energy 
conservation  only  in  average.  More  specifically  the  two 
trajectories  theory88  generated  by  propagating  two  Gaussains,  one 
for  each  electronic  state,  gives  an  ionic  trajectory  and  a 
neutral  one.  If  the  incident  trajectory  is  neutral  then  the 
kinetic  energy  of  the  outgoing  neutral  trajectory  is  the  same  as 
the  kinetic  energy  of  the  incoming  one.  However  the  kinetic 
energy  of  the  outgoing  ionic  trajectory  is  less  than  the  incident 
kinetic  energy  of  the  neutral  by  an  amount  equal  to  the 
ionization  energy.  The  mean  trajectory  method  gives  the  same 
kinetic  energy  for  the  ion  and  neutral.  Furthermore,  the  two 
trajectory  method  can  generate  an  ionic  trajectory  that  is 
trapped  in  the  ionic  well  and  a  neutral  one  which  leaves  the 
surface.  The  mean  trajectory  cannot  split  in  this  way.  It  is 
therefore  not  useful  in  dealing  with  trapping.  Finally,  the  two 
trajectory  method  gives  ion  and  neutral  trajectories  having 
turning  points  on  the  ionic  and  the  neutral  curves,  respectively; 
the  mean  trajectory  has  a  turning  point  at  a  mean  position 
between  the  two  turning  points  mentioned  above.  This  comparison 
between  the  two  methods  shows  the  one  should  not  use  the  MTA  at 
very  low  energy  when  sticking  (i.e.  trajectory  splitting)  is 
important,  or  when  one  intends  to  analyze  detailed  kinetic  energy 
measurements.  However,  it  is  hard  to  believe  that  the  ionization 
probabilities  at  moderate  and  low  energies  are  very  sensitive  to 
the  details  of  the  nuclear  motion.  For  such  situations  the  mean 


trajectory  method  should  be  ar.  adequate  tool  for  calculating  t 
ionization  yield. 
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APPENDIX  A 


The  derivation  of  Eq  (7)  is  similar  to  the  procedure  used 
:c  obtain  the  coherent  state  representation  of  the  path  integral 


-4  H  ( t  .  )  A  t 

<4>  j  I  e  n  J  I  ^ 


<m,.|H(  t  )  ju»  > 
l— At  — — r-2 — - 


-  <^)  j  !  ^  j-l>  exp(“K  <<J,j  I  H ( t )  I  «i» ^ > A t) 


( A- 1 ) 


where  4>--4>(tj).  Here  we  allowed  for  an  explicit  time  dependence 
of  H  and  have  made  the  substitution 


because  this  quantity  already  multiplies  the  small  At.  Also 


=  <4i  j  |<l>  j-<l»  jAt>  =  !4>j>At 


(A-2) 


-  e 


-<4>.j  I  4<  j>At 


Combining  Eqs  (A-l)  and  (A-2)  and  using 


<  41 ,  i  H  j  4i  .  >  =  <  (  Z  x  (tj  |  n>  )  |  H  |  ( Z  x  ( t .. )  !  n>  )  >  with 
J  n  J  n  * 


cf  Eq  (1)  yield  Eqs  (7) -(9). 
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Figure  Captions 


Figure  1 .  The  energies  of  the  neutral  and  the  ionic  states 

involved  in  the ' description  of  the  ionization 

process.  Various  ionic  curves  differ  through  the 

final  energy  e  of  the  transferred  electron,  with 
0£ 

resoect  to  the  Fermi  level .  I  is  the  ionization 
*  00 

potential  of  the  neutral  in  the  absence  of  the 
surface.  6^  is  the  work  required  to  remove  an 
electron  from  the  metal  to  infinity. 

Figure  2.  The  energies  of  the  neutral  and  the  ionic  curves 

invoked  in  discussing  the  neutralization  of  the  io 

The  neutral  curves  all  differ  through  the  initial 

energv  ri  of  the  thermallv  excited  electron  used  f 
P 

neutralizing  the  ion.  & .  is  the  binding  energy  of 


he  ion  to  the  surface. 
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